Excitation lines and the breakdown of Stokes— Einstein relations in supercooled liquids 
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By applying the concept of dynamical facilitation and analyzing the excitation lines that result 
from this facilitation, we investigate the origin of decoupling of transport coefficients in supercooled 
liquids. We illustrate our approach with two classes of models. One depicts diffusion in a strong 
glass former, and the other in a fragile glass former. At low temperatures, both models exhibit 
violation of the Stokes-Einstein relation, D ~ r _1 , where D is the self diffusion constant and r 
is the structural relaxation time. In the strong case, the violation is sensitive to dimensionality 

2/3 n Q 5 

d, going as D ~ r ' for d = 1, and as D ~ r for d = 3. In the fragile case, however, 
we argue that dimensionality dependence is weak, and show that for d = 1, D ~ t~ 0,73 . This 
scaling for the fragile case compares favorably with the results of a recent experimental study for a 
three-dimensional fragile glass former. 



PACS numbers: 64.60.Cn, 47.20.Bp, 47.54.+r, 05.45.-a 



I. INTRODUCTION 

Normal liquids exhibit homogeneous behavior in their 
dynamical properties over length scales larger than the 
correlation length of density fluctuations. For exam- 
ple, the Stokes-Einstein relation that relates the self- 
diffusion constant D, viscosity r/, and temperature T, 



T 

D oc -, 
rj 



(1) 



is usually accurate 0, This relation is essentially 
a mean field result for the effects of a viscous environ- 
ment on a tagged particle. In recent experimental stud- 
ies, it has been reported that the Stokes-Einstein rela- 
tion breaks down as the glass transition is approached 
in supercooled liquid systems 0, 0, IE IE ■ Trans- 
lational diffusion shows an enhancement by orders of 
magnitude from what would be expected from Eq. 
yj, LlU OjjLUJ 113 ■ Here, we show that this breakdown is 
due to fluctuation dominance in the dynamics of low tem- 
perature glass formers. These pertinent fluctuations are 
dynamic heterogeneities [SlIEEHEIIllIllIll. 
Thus, the Stokes-Einstein breakdown is one further ex- 
ample of the intrinsic role of dynamic heterogeneity in 
structural glass formers [22, 123, |24j . 

In the treatment we apply, dynamic heterogeneity is a 
manifestation of excitation lines in space-time |23l | . This 
picture leads to the prediction of dynamic scaling in su- 
percooled liquids, r(l) ~ l z . Here, r(Z) is the structural 
relaxation time for processes occurring at length scale /, 
and z is a dynamic exponent for which specific results 
have been established |23l l24l |25| . This picture and its 
predicted scaling results differ markedly from those de- 
rived with the view that glass formation is a static or 
thermodynamic phenomenon |26l 1271 l28l l29l |3fJ, |3lJ, |32j, 
l33| . It also differs from mode coupling theory which pre- 
dicts singular behavior at non-zero temperature [34L l35j | . 



This paper is organized as follows. In Sec. [H] we in- 
troduce our model for a supercooled liquid with a probe 
molecule immersed in the liquid. Simulation results are 
given in Sees. Hill and II VI Section ITVI also provides ana- 
lytical analysis of the diffusion coefficient and the Stokes- 
Einstein violation, and explains the origin of the decou- 
pling of transport coefficients based on the excitation line 
picture of trajectory space. Comparison of our theory 
with recent experimental results is carried out in Sec.lVl 
We conclude in Sec. IVH with a Discussion. 



II. MODELS 

We imagine coarse graining a real molecular liquid over 
a microscopic time scale (e.g., larger than the molecu- 
lar vibrational time scale), and also over a microscopic 
length scale (e.g., larger than the equilibrium corre- 
lation length). In its simplest form, we assume this 
coarse graining lea ds to a kinetically constrained model 
23. 24, 36, 37, 381 with the dimensionless Hamiltonian, 



AT 



H = '^2n l , (m = 0, 1). 



(2) 



Here, n; = 1 coincides with lattice site i being a spatially 
unjammed region, while rii = coincides with it being 
a jammed region. We call rii the "mobility field" . The 
number of sites, N , specifies the size of the system. From 
Eq. (|5J, thermodynamics is trivial, and the equilibrium 
concentration of defects or excitations is 



(rii) 



1 



1 + exp(l/T) 



(3) 



where T is a reduced temperature. We make explicit 
connection of T with absolute temperature later when 
comparing our theory with experimental results. 
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The dynamics of these models obey detailed balance 
and local dynamical rules. Namely, 

rii = < rii = 1, (4) 

where the rate constants for site i, fcj and k\ \ depend 
on the configurations of nearest neighbors. For example, 
in dimension d = 1, 

= e- 1 ' f f(jn.. 1 ,n i+1 ), (5) 
= f(ni-i,n i+ i), (6) 

where /(ni_i,rij+i) reflects the type of dynamical facil- 
itation. In the Frcdrickson-Andersen (FA) model 36], a 
state change is allowed when it is next to at least one 
defect. The facilitation function in this case is given by, 

/faO^-i, th+i) = rii-i + n l+ i - m-iUi+i. (7) 

In the East model 37], dynamical facilitation has direc- 
tional persistence. The facilitation function in this case 
is 

/Eafrt(ni_l,n i+1 ) = m-i. (8) 

In order to study translational diffusion in supercooled 
liquids, we extend the concept of dynamic facilitation to 
include a probe molecule. The dynamics of a probe will 
depend on the local state of the background liquid. When 
and where there is no mobility, the diffusive motion of 
the probe will be hindered. When and where there is 
mobility, the probe molecule will undergo diffusion easily. 
As such, in a coarse-grained picture, the probe molecule 
is allowed to jump from lattice site i to a nearest neighbor 
site when site i coincides with a mobile region, n% = 1. In 
order to satisfy detailed balance, we further assume that 
the probe molecule can move only to a mobile region, i.e., 

x(t + 8t) = x(t) ± 8x ■ n x ■ n x ±sx, (9) 

where x(t) denotes the position of the probe at time t. 
Units of time and length scales are set equal to a Monte 
Carlo sweep and a lattice spacing, respectively. 

III. COMPUTER SIMULATIONS 

Using the rules described in Sec.[nJ we have performed 
Monte Carlo simulations of diffusion of a probe molecule 
in the FA and East models for various temperatures. For 
the purpose of numerical efficiency, we have used the con- 
tinuous time Monte Carlo algorithm [3{J |4jj . In the all 
systems, N was chosen as N = 100/c, and the simula- 
tions were performed for total times T w lOOr, with r 
being the relaxation time of the model. Averages were 
performed over 10 3 to 10 5 independent trajectories. 
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FIG. 1: Typical trajectories of a probe molecule in one- 
dimensional models. The probe molecule (black line) under- 
goes a diffusive process in the trajectory space that consists 
of gray (mobile) and white (immobile) regions, (a) FA model 
at f = 3; (b) FA model at f = 0.8, and (c) East model at 
T = 0.8. 



In Fig. ^ we show typical trajectories of probe 
molecules in the FA and East models. In the high temper- 
ature case, trajectory space is dense with mobile regions 
and there are no significant patterns in space-time. As 
such, the dynamics is mean-field like. It is for this reason 
that the relaxation time in this case is inversely propor- 
tional to the equilibrium probability of excitation, c (see, 
for example, Ref. 01 )• The probe molecule executes dif- 
fusive motion, without being trapped in immobile regions 
for any significant period of time. 

The low temperature dynamics is different. Mobility 
is sparse, defects tend to be spatially isolated at a given 
time, and trajectory space exhibits space-time patterns. 
See Fig. db) and (c). Because of the facilitation con- 
straint, an immobile region needs a nearest mobile region 
to become mobile at a later time. The excitations there- 
fore form continuous lines and bubble-like structures 
in trajectory space. While inside a bubble, the probe 
molecule will be immobilized. See, for example, the seg- 
ment of the trajectory of a probe molecule for < t < 500 
in Fig. ^b). Due to exchanges between mobile and im- 




FIG. 2: Mean squared displacements of the probe molecules 
are shown for the three different cases illustraed in Fig. Q 




FIG. 3: Diffusion coefficients for the FA and East models as 
functions of 1/T. Dashed lines are a guide to the eye. 



mobile regions, an immobile region can become mobile 
after a period of time. At that stage the probe molecule 
can perform a random walk until it is again in an immo- 
bile region. The motion of a probe molecule will manifest 
diffusive behavior over a time long enough for many dy- 
namical exchanges to occur. In the East model at low 
temperatures such as pictured in Fig. Uc), the bubbles 
in space-time form hierarchical structures |23|. 

Figure [3 plots mean square displacements of probe 
molecules for the FA and East models for three differ- 
ent cases pictured in Fig. ^ In the high temperature 
case, the mean square displacement reaches its diffusive 
linear regime after a very short transient time. In the low 
temperature case, the probe molecule in the East model 
case reaches the diffusive regime after a longer time and 
over a larger length scale than that in the FA model with 
the same reduced temperature. 

IV. STOKES— EINSTEIN VIOLATION 

A. Diffusion Coefficient 

Figure |3| plots the diffusion coefficient of a probe 
molecule for the FA and the East models. The diffusion 
coefficient is determined from the mean square displace- 
ment, 

D = ta M>, do) 

t^oo t 

where Ax(t) = x(t) — x(0). Error estimates for our sim- 
ulations are no larger than the size of the symbols. 

In the FA model, the diffusion coefficient exhibits Ar- 
rhenius behavior for T < 1. This behavior reflects the 
fact that relaxation dynamics in the FA model is similar 
to that of a strong liquid. In this regime, over more than 



4 orders of magnitude in D, the slope of \ogD vs T _1 
is close to 2. This result is consistent with the expected 
low temperature scaling, 

D FA ~ c 2 ~ exp(-2/f ), (11) 

as discussed in the next subsection. In the East model 
case, also pictured in Fig. |3J the diffusion coefficient 
decreases more quickly than Arrhenius. This super- 
Arrhenius behavior is due to the hierarchical nature of 
dynamics in the East model [42|. 

Comparing the diffusion coefficients with the relax- 
ation times of the background liquids demonstrates 
Stokes-Einstein violation in both models. The relax- 
ation times, t, of the FA and the East models at dif- 
ferent temperatures have been determined in prior work 
[43L liij. When the Stokes-Einstein relation is satisfied, 
Dt ~ const. This behavior occurs in the FA and East 
models when T > 2, but Fig. 21 shows that Dt is en- 
hanced from that behavior by 2 or 3 orders of magnitude 
when T < 1. Bear in mind, these deviations from Stokes- 
Einstein are d = 1 results. The appropriate generaliza- 
tion of the FA model to d — 3 does not exhibit such large 
deviations. On the other hand, we expect that gener- 
alizations of the East model, which is hierarchical and 
therefore fragile, will have weak dimensional dependence 
and continue to exhibit large deviations for d = 3. We 
turn to the arguments that explain these claims now. 

B. Scaling Analysis 

For high temperatures, the local mobility field will tend 
to be close to its mean value, c. As such, both the re- 
laxation mechanism of the material and the diffusional 
motion of the probe molecule make use of the same local 
mobility fields. For this reason, the diffusion coefficient 
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FIG. 4: Violations of the Stokes-Einstein relation are similar 
in the d — 1 FA and East models. Dashed lines are a guide 
to the eye. 



and the relaxation time scale are strongly coupled in this 
regime, leading to the Stokes-Einstein relation. 

At low temperatures, however, the dynamics of the 
system is not so simply related to the mean mobility 
field. Here, the fluctuations of bubble structures dom- 
inate. The relaxation time of the background liquid will 
approximately scale as the longest temporal extension of 
bubbles. The persistence time of an individual lattice 
site, ipersj is the time for which that site makes its first 
change in state. Its typical size will be intimately tied to 
the structural relaxation time of the liquid. For the FA 
model in d = 1 , 



T ~ (tpcrs) ~ C- 3 . (12) 

See, for example, Refs. |2j ^ . 

This result is consistent with a simple argument con- 
cerning diffusive motions of excitation lines in the low 
temperature FA model 0]. In particular, the struc- 
tural relaxation times in the FA model is given by the 
time in which a typical bubble structure looses its iden- 
tity through wandering motions of excitation lines. The 
excitation line has a local diffusivity of T> ~ c. (We 
use caligraphic T> to distinguish this diffusion constant 
for excitations from that for particles, D.) In order to 
form a bubble, an excitation line needs to wander dis- 
tance of the order of the typical length between defects, 
^cq ~ c _1 . Therefore, the mean relaxation time is given 
by r ~ IIJV ~ a" 3 . 

When the probe molecule is at the boundary of a bub- 
ble, it may not need to wait until the bubble closes in 
order to undergo diffusion; rather, it can remain within 
mobile cells and diffuse around the boundary of the bub- 
bles. In this way, translational diffusion will be more 
facilitated than structural relaxation, leading to an en- 
hanced diffusion in the fluctuation dominated low tem- 
perature region. Specifically, consider the dynamical ex- 




FIG. 5: A section of Fig. ^b) illustrating the meanings of 
exchange times, to and ti. to is a time a site spends in a 
bubble, and t\ is a time it spends in a surrounding boundary. 



change times, i.e., 
a given lattice site, 
for an ?ii = state 
n.j = 1 state. The 
in a mobile region 
ment of the probe 
diffusive steps that 
trajectory, N, 



the times between flipping events for 
See Fig.|3J to is such a time duration 
and t\ is such a time duration for an 
probe molecule can move only while 
. Further, the mean square displace- 
will be proportional to the number of 
a probe molecule will take during the 



([Ax(t)] 2 )~Af. 



T 



(to) + (h) 



(13) 



Here, T is the length of a long trajectory in the FA model. 
The average duration of the defect state, (ti), is inversely 
proportional to the probability of a lattice site being mo- 
bile, c, times the flip rate, k\~\ Since k^' ~ 0(1), we 
have, 



<*!> 



From detailed balance, therefore, 



(to) 



(14) 



(15) 



Since (t\) <C (to) in the low temperature region, 
Eqs. (JT3J-0 give 



((Ax) 2 ) 1 2 
d F a ~ y c 



(16) 



This result explains Eq. {HJ . Together with Eq. (O , it 
leads to 



D 



FA 



(17) 



with £ = 2/3 in the d = 1 FA model case. This scaling is 
to be contrasted with the Stokes-Einstein result, f = 1. 

Numerical simulation [24| and renormalization group 
analysis |2J| of higher dimension generalizations of the 
FA model indicate that for d = 3, r ~ c~ 21 . However, 
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FIG. 6: Scaling of Stokes-Einstein violation in d = 1. Circles 
and squares indicate computed results for the FA and East 
models, respectively. 



the scaling D ~ c 2 remains true for all dimensions as 
it is based solely on detailed balance. Thus, for d = 3, 
£ ~ 0.95. In other words, there is only a weak breakdown 
in the Stokes-Einstein relation for strong liquids in d = 3. 

In the East model case, both the diffusion coefficient 
and the relaxation time show super- Arrhenius behavior. 
The hierarchical, fractal structure of pattern develop- 
ment in trajectory space for the East model does not 
allow a simple scaling analysis of the diffusion coefficient, 
and it is not obvious whether temperature independent 
scaling exists. One can define a temperature dependent 
scaling exponents, a(T) and z(T), 



so that 



D 



^East 
T 

East ^ 



c a( - T \ 



-a(T)/ Z (T) 



(18) 
(19) 

(20) 



Interestingly, our numerical results indicate that £ = 
a/z rs 0.73 is independent of temperature as shown in 
Fig-El This exponent, £ w 0.73 for the d = 1 East model, 
is very close to what many experiments and simulations 
have found for three-dimensional glass forming liquids. 
For example, a recent experiment finds that £ « 0.77 in 
the self-diffusion of iris-naphthylbenzene(TNB) It 
was found that £ w 0.75 in a molecular dynamics sim- 
ulation of Lennard-Jones binary mixture 45] and a re- 
cent detailed scaling analysis of numerical results shows 

£ « 0.65 m. 

Presumably, such good agreement of scaling relation 
between the d = 1 East model and higher dimension sys- 
tems arises due to directional persistence of facilitation in 
the fragile liquid [23L |24| . This persistence in higher di- 
mensions causes motion to be effectively one-dimensional 



|24|. Therefore, dimensionality is not very significant for 
fragile glass formers. As such, for fragile systems, we 
expect that the scaling relation of the Stokes-Einstein 
violation will be reasonably well described by the d = 1 
East model. Based on this expectation, we further purse 
the comparison between theory and experiment. 



V. COMPARISON WITH EXPERIMENT 

Swallen et al. || measured the self-translational diffu- 
sion coefficient of TNB near the glass transition temper- 
ature. They observed an increase of Drj/T from its high 
temperature limit by a factor of 400 near the glass tran- 
sition temperature. In order to compare our results with 
these experiments, we need to determine the excitation 
concentration, c, as function of temperature. Since TNB 
behaves as a fragile liquid, we determine the excitation 
concentration as a function of temperature by fitting the 
viscosity data of TNB (4^ with the generalization of the 
East model formulas to higher dimensions [24j . Namely, 



lnr 



1 



rfln2 



MsAOf , 



(21) 



where g is the number of equally likely persistence direc- 
tions on a cubic lattice, and 



ln(c) = ln( Cfl ) - J 



1 



1 



T Tn 



(22) 



The parameter J is the energy scale associated with cre- 
ating a mobile region from an immobile region, and Tr 
is an appropriate reference t emp erature. Details on the 
fitting can be found in Ref. 24]. Taking g = 8 (the cu- 
bic lattice value) and Tn as the temperature at which 
logr is half the value of log r(T g ), we determine that 
J/T g w 21.7, and \og w {c R /g) w -1.28. The reduced 
temperature, T of the East model is related to absolute 
temperature by 



(23) 



Once we have determined the excitation concentration 
as a function of the temperature, we can compare exper- 
imental data with our computed results for the Stokes- 
Einstein violation in the East model case. Based on the 
argument that the scaling relation of the Stokes-Einstein 
violation (D ~ r~ remains robust in higher dimensions 
and from the dimensional dependence of Eq. (|21|l . we 
expect 



m(£>T)U =3 w i In ODr) | d=1 . 



(24) 



In Fig. we use this relationship to compare the extent 
of the Stokes-Einstein violation of the experimental sys- 
tem with our East model results. The agreement between 
the two appears excellent. 
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FIG. 7: Comparison between the East model prediction and 
experiments on supercooled TNB, Ref. 



VI. DISCUSSION 

There have been previous theoretical studies on the 
violation of the Stokes-Einstein relation in supercooled 
liquid systems. For example, Kivelson and Tarjus have 
argued that the Stokes-Einstein violation can be under- 
stood from their "frustration-limited domain" model for 
supercooled liquids 0, . Assuming a distribution of 
local relaxation times associated with domain structures, 
this model describes the translational diffusion and vis- 
cosity as corresponding to different averaging process of 
such a distribution. Their idea contrasts to ours in that 
the domain structure in their work is purely static, and 
the exchange between different domains are not consid- 
ered. 

Hodgdon and Stillinger have proposed a fluidized do- 
main model |sL llO^ . In their work, it is assumed that the 
system consists of a sparse collection of fluid-like domains 
in a background of more viscous media, and fluid-like do- 
mains appear and disappear with a finite life-time and 



rate. Relaxation times are determined by the rate of ap- 
pearance of the fluid-like domains, while translation dif- 
fusion also depends on the life-time of the domains. To 
the extent that these domains refer to space-time and 
not simply space, this picture is not inconsistent with 
ours. Xia and Wolynes have applied the so-called "ran- 
dom first order transition theory" |3^| to the Hodgdon- 
Stillinger model I n this case, the picture is both 

mean field and static and decidedly contrary to our fluc- 
tuation dominated and dynamic view. 

From the perspective that Stokes-Einstein violation is 
a manifestation of fluctuation dominated dynamics, one 
expects that similar decoupling behavior occurs between 
other kinds of transport properties near the glass tran- 
sition. The extent to which such decoupling can appear 
depends upon microscopic details in the specific trans- 
port properties and materials under study. For exam- 
ple, molecular rotations of a probe will be coupled to the 
mobility field, but less so than translations. Indeed, sin- 
gle molecule experiments indicate that rotations persist 
in both mobile and immobile regions of a glass former 
[IM l49l IHoj . Rotational motions can therefore average 
the effects of dynamic heterogeneity to a greater extent 
than translational motions. As such, decoupling of rota- 
tional relaxation from structural relaxation can be more 
difficult to detect than violations of the Stokes-Einstein 
relation. Precisely how such effects might be detected 
seems worthy of further theoretical analysis. 
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